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Abstract 

We consider rotational initial data for the two-dimensional incompress- 
ible Euler equations on an annulus. Using the convex integration frame- 
work, we show that there exist infinitely many admissible weak solutions 
(i.e. such with non-increasing energy) for such initial data. As a conse- 
quence, on bounded domains there exist admissible weak solutions which 
are not dissipative in the sense of P.-L. Lions, as opposed to the case with- 
out physical boundaries. Moreover we show that admissible solutions are 
dissipative provided they are Holder continuous near the boundary of the 
domain. 

Dedicated to the memory of Professor Mark Iosifovich Vishik 



1 Introduction 

The study of weak solutions of the incompressible Euler equations is motivated 
by (at least) two aspects of fluid flow: the presence of instabilities, most notably 
the Kelvin-Hclmholtz instability, and fully developed 3-dimensional turbulence. 
Concerning the latter, an important problem arises in connection with the fa- 
mous 5/3 law of Obukhov-Kolmogorov and the conjecture of Onsager regarding 
energy conservation. We refer to [3j[15] and [6l[9j[19] for more information and 
recent progress regarding this problem. 

Concerning the former, it has been the subject of intensive research to de- 
fine a physically meaningful notion of weak solution, that can capture the basic 
features of such instabilities and be analytically well behaved at the same time. 
Due to the lack of an analogous theorem to the existence of Leray-Hopf weak 
solutions of the Navier-Stokes equations, several weaker notions have been con- 
sidered. 

Dissipative solutions of the incompressible Euler equations were introduced 
by P.-L. Lions [52] as a concept of solution with two desirable properties: (i) 
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existence for arbitrary initial data, and (ii) weak-strong uniqueness, meaning 
that a dissipative weak solution agrees with the strong solution as long as the 
latter exists. Dissipative solutions have been shown to arise, among others, as 
viscosity [22] or hydrodynamic [H] limits of the incompressible Eulcr equations. 
The major draw-back of dissipative solutions is that, in general, the velocity 
field does not solve the Eulcr equations in the sense of distributions. 

Weak solutions (i.e. distributional solutions with some additional properties) 
on the other hand have been constructed by various techniques, see [TUIITT1|2"51 - 
H71HII1IH3I33] ■ Many of these results come with a high level of non- uniqueness, 
even violating the weak-strong uniqueness property - we refer to the survey [12] . 
In particular, in |33j the existence of global in time weak solutions was shown 
for arbitrary initial data. 

Due to the high level of non-uniqueness, a natural question is whether there 
are any selection criteria among weak solutions. With this regard, it has been 
noted in [TUT^] that, in the absense of boundaries a weak solution is dissipative 
in the sense of Lions, provided the weak energy inequality 



holds. In [TT] this condition is referred to as an admissibility condition, in 
analogy with the entropy condition used in hyperbolic conservation laws [8|. 
Admissibility turned out to be a useful selection criterion among weak solutions, 
since already in the weak form in ([!]) it implies the weak-strong uniqueness 
property of dissipative solutions (stronger versions of the energy inequality are 
discussed in [H]). This is even the case not just for distributional solutions but 
also for measure- valued solutions, see [5]. 

Despite the weak-strong uniqueness property, there exists a large, in fact L 2 
dense set of initial data on the whole space or with periodic boundary condi- 
tions |31j (see also [29]), for which the initial value problem admits infinitely 
many admissible weak solutions. Such initial data, called "wild initial data" , 
necessarily has to be irregular. 

The non-uniqueness of admissible weak solutions is intimately related to the 
presence of instabilities. For instance, in |30j the non-uniqueness of admissible 
weak solutions was shown for the flat vortex sheet initial data 



extended periodically to the torus T d . Note that the stationary vector field is 
an obvious solution in this case, but the statement in [30| is that there exist 
infinitely many non-stationary solutions. A common feature in these solutions 
is that for time t > they exhibit an expanding "turbulent" region around the 
initial vortex sheet, much akin to the propagation of singularity in the classical 
Kclvin-Hclmholtz problem. Further examples of this nature appeared in [4] and 
recently in [7] for the compressible Euler system. 




for almost every t > 
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Motivated by the idea that it is the underlying Kelvin-Helmholtz instability 
that is responsible for the non-uniqueness of admissible weak solutions, we study 
in this note the case of domains with boundary. We show that the presence of 
a (smooth) boundary can lead to the same effect of an expanding turbulent 
region as in |30j . As a corollary, we observe that admissibility docs not imply 
the weak-strong uniqueness property in domains with boundary. 

2 Statement of the main results 
2.1 Formulation of the equations 

We study weak solutions of the initial and boundary value problem for the 
incompressible Euler equations 

d t v + v ■ W + Vp = 

divw = (3) 
v\ t =o = v 

complemented with the usual kinematic boundary condition 

v\dn - v = 0. 

Here, Q c W l , d > 2, is a domain with sufficiently smooth boundary, T > a 
finite time, v : Q x [0, T) -> R d the velocity field, p:f!x(0,T)->l the scalar 
pressure, vq the initial velocity and v the inner unit normal to the boundary of 
0. 

In order to give the precise definition of weak solutions, consider the space 
of solenoidal vectorfields on S7 (cf. Chapter III of |T5] ) , 

ff(fi) ={veL 2 (n;R d ) : [ v-Vpdx = 
Jn 

for every p e W^(fl) such that Vp e L 2 {ti)}. 

Let i>o 6 H(Cl). An admissible weak solution of © with initial data vq is 
defined to be a vectorfield v € L°°(0,T;H(S1,)) such that for every test function 
4> e C c °°(fi x [0,T);R 2 ) with div0 = 0, we have 

T r 

(dt<t> • v + V '<j> ■ v ® v) dxdt + / vq(x) ■ <j>{x, 0)dx = 0, 



Jn 

and the energy inequality (flj holds. 

We remark in passing that in fact one may assume that admissible weak 
solutions are in the space C([0, T); H w (il)), where H w (il) is the space if(fi) 
equipped with the weak L 2 -topology. Indeed, dissipativc solutions of Lions are 
also defined in this space. Nevertheless, for simplicity we will just treat the 
velocity fields as elements in the larger space L°°(0,T; iT(Q)). 
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2.2 Rotationally symmetric data 

In the present paper, we consider rotationally symmetric initial data in two 
dimensions. It should be noted that the restriction to 2 dimensions is purely 
for simplicity of presentation - the constructions and the methods can be easily 
extended to higher dimensions. Similarly, we will consider as domain an annulus 
purely for simplicity of presentation - the nontrivial topology of the domain does 
not play a role in our results. 

By "rotational" we mean initial data of the form 

vq(x) = ao(r)(sm9, -cos#) (4) 

on an annulus 

n = {xeR 2 :p<\x\<R}, (5) 

where < p < R < oo. Vector fields as in (0| are known to define stationary 
solutions to the Euler equations regardless of the choice of ao, and are frequently 
used as explicit examples in the study of incompressible flows p~lf23 . 28 . 

Fix a radius ro with p < tq < R and consider the initial data on the annulus 
given by ^ with 

Mr) -If 'ZT^R (6) 

which corresponds to a rotational flow with a jump discontinuity on the circle 
{r = r }. 

Theorem 1. Let Q be an annulus as in {5P, T > a finite time, and vq be 
rotational as in ^ and ^\). Apart from the stationary solution v(-,t) = vq, 
there exist infinitely many non- stationary admissible weak solutions of the Euler 
equations on O x (0,T) with initial data vq. Among these, infinitely many have 
strictly decreasing energy, and infinitely many conserve the energy. 

Our proof, given in Section 2] below, relies on the techniques from QT| and 
is similar to the construction in [30] . 

Regarding the quest for suitable selection principles, a much-discussed crite- 
rion is the viscosity solution, defined to be a solution obtained as a weak limit 
of Leray-Hopf solutions as viscosity converges to zero. In the case of the initial 
data in ([2]) it is an easy exercise (see for instance [4]) to show that the viscosity 
solution agrees with the stationary solution. In the rotational case ([6]) the same 
is true, as wc show in Section [5] below: 

Proposition 2. Let fi c R 2 be an annulus and let initial data be given by Q). 
Then every sequence of Leray-Hopf solutions of the Navier-Stokes equations with 
viscosities tending to zero which correspond to this initial data will converge 
strongly to the stationary solution v(-,t) = Vq of the Euler equations. 

Finally, we discuss the relation between admissible weak solutions and dis- 
sipative solutions of Lions in bounded domains. For the convenience of the 
reader we recall in Section |6] the precise definition of dissipative solutions. As a 
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corollary to Theorem [T] we show in Section [5] that, contrary to the case without 
boundaries, admissible weak solutions need not be dissipative: 

Corollary 3. On Q there exist admissible weak solutions which are not dissi- 
pative solutions. 

Corollary|3]says that in the presence of boundary the weak-strong uniqueness 
might fail for admissible weak solutions. On the technical level the explanation 
for this lies in the observation that the notion of strong solution in a bounded 
domain does not allow any control of the boundary behaviour. Therefore in 
Scction[7]wc study what happens when additional boundary control is available: 

Theorem 4. Let f2 c M. 2 be a bounded domain with C 2 boundary. Suppose v is 
an admissible weak solution of on Q for which there exists some S > and 
a > such that v is Holder continuous with exponent a on the set 

T s = {xeTi:dist(x,dri) < 5} , 

uniformly in t. Then v is a dissipative solution. 



3 Subsolutions and convex integration 

In order to prove Theorem [T] we recall the basic framework developed in pTJlfTT] . 
with slight modifications to accomodate for domains with boundary. For further 
details we refer to the survey [H] and the recent lecture notes |29| . 

To start with, recall the definition of subsolution. To this end let us fix a 
non-negative function 

e€L oo (0,T;L 1 (O)), 

which will play the role of the (kinetic) energy density. We will work in the 
space-time domain 

n T ■= n x (0,T), 

where O c R d is either an open domain with Lipschitz boundary or = T d . 

Definition 5 (Subsolution). A subsolution to the incompressible Eulcr equa- 
tions with respect to the kinetic energy density e is a triple 

with v e L°°(0,T;H(Q,)), u e L\ JQ, T ), q e £»'(f2 T ), such that 

d t v + divw + Vo = . c t i *i j* / n \ 

,. _ _ m the sense of distributions; (7 

divv = 0, 

and moreover 

u®w-u<|e/ a.e. (x,t). (8) 
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Here Sq denotes the set of symmetric traceless d x d matrices and I is 
the identity matrix. Observe that subsolutions automatically satisfy ^\v\ 2 < e 
a.e. If in addition © is an equality a.e. then v is a weak solution of the Eulcr 
equations. 

A convenient way to express the inequality is obtained by introducing 
the generalized energy density 

, N d, _ , 

e{v,u) = -\v ®v - uloo, 

where | • |oo is the operator norm of the matrix (= the largest eigenvalue for 
symmetric matrices). The inequality ([8]) can then be equivalently written as 

e(v,u) < e a.e. (9) 

The key point of convex integration is that a strict inequality instead of ([8]) 
gives enough room so that high-frequency oscillations can be "added" on top 
of the subsolution - of course in a highly non-unique way - so that one obtains 
weak solutions. It is important also to note that, since in the process of convex 
integration only compactly supported (in space-time) perturbations are added 
to the subsolution, the boundary and initial conditions of the weak solutions 
so obtained agree with the corresponding data of the subsolution. This is the 
content of the following theorem, which is essentially Proposition 2 from . 

Theorem 6 (Subsolution criterion). Let e e L°°(Qt) and (v,u,q) be a sub- 
solution. Furthermore, let U c fl T a subdomain such that (v,u,q) and e are 
continuous on U and 

e(v,u) <e on U 

_ n (io) 

e(y,u) = e a.e. VLt^U 

Then there exist infinitely many weak solutions v e L°°(0,T;iJ(£l)) of the Euler 
equations such that 

v = v a.e. Qt n U, 
\\v^ = e a.e. Qt, 
p = q--;e. a.e. J7t- 

// in addition 

v(-,t) - v (-) in L 2 (n) as t ^ 0, (11) 
then v solves the Cauchy problem ([3]). 

We also refer to [53] , where a detailed discussion of the convex integration 
technique can be found - in particular the above theorem is Theorem 7 of |29j . 
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4 Non-Uniqueness for Rotational Initial Data 

In this section we wish to apply the framework of Section [3] to prove Theorem 
[TJ Thus, we set 

!l:={i;ei 2 : p < \x\ < R} 
to be an annulus, fix ro € (p, R) and set 

f-rWa; 1 Id < rn, 
N3 o, 
j^x 1 \x\>r , 

where x 1 = | X2 I . We will construct subsolutions by a similar method as in [30] • 

V"* 1 / ... 
Owing to Theorem [S] of the previous section, it suffices to show the exis- 
tence of certain subsolutions. We fix two small constants A > ("turbulent 
propagation speed") and e > ("energy dissipation rate"), to be determined 
later. 

We look for subsolutions (v,u,q) (c.f. Definition [5] - the energy density 
function e is still to be fixed) of the form 

where a(r, 0) = ao(r) and (r, 9) denotes polar coordinates on R 2 , 

U ^ X ' >~ \ sm9 -cosd )[ 7 (r,t) -0(r,t) ){ sind -cos9 J 



_ I /3cos(20)+7sin(20) /3sin(26») - 7 cos(26») \ 
~\f3sm(29)- 1 cos{29) -/?cos(26>) - 7 sin(2(9) J 



(13) 



and 

q = q(r). 

As a side remark, note that the choice a(r, t) = ao(r) for alH > 0, (3 = -\c? , 
7 = 0, and 

q{r) = \c? + f^ds (14) 

yields the well-known stationary solution (the integral in the formula for q rep- 
resents the physical pressure). 

We insert this ansatz into ([7]) to arrive at two equations. More precisely, 

,! n i /cos#\ , „ i /-sin#\ , 

using tnc formulas V x r = . „ and V x " = - n , we obtain 
6 y$,m9 J r \ cos9 J' 



d t a sin 9 + d r (3 [cos 9 cos(26>) + sin 9 sin(2(9)] + d r ~f [cos 9 sin(26>) - sin 9 cos(2#)] 
2 

+ -j3 [sin 9 sin 
r 

+ d r qcos9 = 



2 2 
+ -/3 [sin 9 sin(26») + cos (9cos(2(9)] + -7 [-sin6>cos(26>) + cos6>sin(2#)] 
r r 
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and 

-d t acos6 + d r P [cos 6 sin(26>) - sin0cos(20)] + <9 r 7 [-008 6*008(20) - sin0sin(20)] 

+ -j(9[-sin0cos(20) + cos0sin(20)] + -7 [-sin0sin(20) - cos0cos(20)] 
r r 

+ d r qsm8 = 0. 

If we multiply the first equation by sin 9 and add it to the second one multiplied 
by cos0, use the identities cos 2 0-sin 2 = cos(20) and 2sin0cos0 = sin(20), 
and then separate by terms involving sin(20) and cos(20), respectively, we will 
eventually get the two equations 

d r (3 + -(3 + d r q = 



d t a + d r j + -7 = 



2 ( 15 ) 



It can be easily verified that these equations are equivalent to the original system 
(|7|) for our ansatz. 

If we set q(r) as in (fbfj) and (3 = -\a 2 , the first equation will be satisfied, 
in nice analogy with [30] (up to a sign). Also, the second equation is similar 
to [30], but it involves the additional "centrifugal" term ^7. Therefore, we 
cannot simply set 7 = |a 2 as in [3D] to obtain Burgers' equation. However, 
observing that d r (r 2 ^f) = 2rj + r 2 9 r 7, we set 

1 
r 

and 



a(r,t) = -f(r,t) 



^<W*)-^-,V), (16) 

so that the second equation in (|I5|) . after multiplication by r 2 , turns into Burg- 
ers' equation 

dtf + \d r {f)=0. (17) 



The initial data ([6]) for a then corresponds to 

/(r,0) 



-1 if p < r < ro 
1 if r < r < R. 



Then, for this data, Burgers' equation (|T7)) has a rarefaction wave solution for 
t e [0,T], provided A > is sufficiently small (depending on T and p < ro < R), 
which can be explicitly written as 



f(r,t) 



-1 if p < r < ro - Xt 

if r - Xt < r < r + Xt (18) 
1 if ro + Xt < r < R. 
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Therefore, by setting a(r,t) = ~2f(r,t) for / as in (JTHJ), (3 = -\o? , 7 as in 
(fT6|) , and q as in (| 14[) , we obtain a solution of the equations (O with initial data 
corresponding to (|T^|) . 

It remains to study the generalized energy. Since u is given by (|13p and 
moreover 

, ,,9 / cos 2 - sin 6* cos 9 \ 
v®v = a(r,t) _ sin0cos0 CQS 2 d J 

I cos9 sm9 WO W cos6» sin0 \ 
~\ sm9 -cose )[ a(r,i) 2 J \ sinfl -cosfl J' 

and since the eigenvalues of a matrix are invariant under conjugation by an 
orthogonal transformation, in order to determine e(v,u) = \v®v — u\oo it suffices 
to find the largest eigenvalue of 

1-7 « 2 ^j"U(^- 2 " 2 ) W I' 

It is easily calculated, taking into account \a\ < -j and A > 0, that 

1 9 A / 1 



e(«,«) = -a + 2^" r a J 



-l[l-(l-, 2 A)(l-/(r,;) 2 )] 



(19) 



Finally, we set 



e -(r,t) = ^ I [l- £ (l-r 2 A)(l-/(r,t) 2 )], 

where e is sufficiently small so that e > 0. Observe that 

e(u, u) < e < — | 2 inilx- 

More precisely, we have the following result, summarizing the calculations in 
this section: 

Proposition 7. For any choice of constants e, A satisfying 

1 r -p R-r 



< A < min 
0<e< 



i? 2 ' T 



l-p 2 A 

there exists a subsolution (w, u, q) in J7t with respect to the kinetic energy density 
e(r,t) = ^[l-e(l-r 2 A)(l-/(r,t) 2 )] 
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and with initial data v(x,0) = vq(x) from (|12p . such that, with 
U := {x € R 2 : r Q -Xt< \x\ < r + At} 

we have 



e(v, u) < e m U, 
e{v, u) = e m ily \ £Y. 



We can now conclude with the proof of Theorem [T] 

Proof of Theorem^ We apply Proposition [7] above with e > to obtain a sub- 
solution (v,u,q). According to Theorem [5] with this subsolution, there exist 
infinitely many weak solutions v e L°° (0,T; H(il)) such that |w| 2 = 2e almost 
everywhere in 0^ and with initial data vq. To check that these are admissible, 
observe that 

/ \v(x, t)\ 2 dx = [ 2e(x, t) dx < — — dx = [ \vo(x)\ 2 dx. 
Jn Jn Fr Jn 

Finally, observe that we obtain strictly energy-decreasing solutions by choosing 
e > and energy-conserving solutions for e = 0. □ 

5 Uniqueness of the Viscosity Limit 

Proof of Proposition [H Consider the Navier-Stokes equations with viscosity e > 
0: 

dtv e + v e • X7v e + Vp c = eAv e 
div v e = 
v t {-,0) = v a 
v t \dn = 0. 

It is known that the Navier-Stokes equations in two space dimensions admit 
a unique weak solution (the Leray-Hopf solution) which satisfies the energy 
equality 

- / \v e (x,t)\ 2 dx + e / / \Vv e (x, s)\ 2 dxds = - / \vo(x)\ 2 dx 
2 Jn Jo Jn 2 Jq 

for every t 6 [0,T], see e.g. [17] for details. It turns out that if the initial data 
vq has the rotational symmetry in (j4]), then the (unique) Leray-Hopf solution 
will have the same symmetry. 

To show this, we take the ansatz 

v e (x,t)=a e (r,t)l S[ ^ a ) (21) 
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and p e = p e (r), again using polar coordinates. Insertion of this ansatz into the 
first equation of (|20| yields 



a 2 



dta e sin 9 cos 9 + d r p e cos 9 



' d r a t 2 a e \ 

+ oia f -\sm0 . 



If we choose 



Pe{r) = / as 



and divide by sin (9, we end up with the parabolic equation 



I d r a e 2 
o t a t = e + o r a e 



Insertion of our ansatz into the second equation of (f20f also gives (|22|) . as one can 
easily check by a similar computation. Moreover, the divergence-free condition 
is automatically satisfied, the initial condition becomes 

a £ (-,0)=a o (23) 

with ctQ defined by (|6|), and the boundary condition translates into 

a e (p) = a e (R) = 0. (24) 

Thus we obtain the well-posed parabolic initial and boundary value problem 
([22]l. ([23]l. (|24]l. By well-known results (cf. e.g. [14], Section 7.1), this parabolic 
problem admits, for each e > 0, a unique weak solution. But our calculations so 
far show that, if a e is a solution to the parabolic problem, then the corresponding 
v e defined by (|2Tj) is the (unique) Leray-Hopf solution of the Navier-Stokcs 
problem (|20[) . and at the same time it satisfies the initial and boundary value 
problem for the heat equation: 

d t v e = eAv c 
div v e = 
v e (-,0) = v Q 
v £ \dn = 0. 

Since the solutions of the heat equation converge strongly to the stationary 
solution, and since we have shown that for our particular initial data the heat 
equation coincides with the Navier-Stokcs equations, the proposition is thus 
proved. □ 

Remark 8. The previous discussion can be extended to initial data on a cylinder 
of the form Z = £!xTcR 2 xT, where ft c M 2 is still the annulus. Indeed, for so- 
called 2 1/2 dimensional initial data Vq(xi,X2) = (vo(xi,X2),w(xi,X2)) on Z, 
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where "Uo is as in there may exist infinitely many admissible weak solutions, 
but only the solution given by 

V(xi,x 2 ,t) = (v (xi,x 2 ),w(xi - (v )it,x 2 - (v ) 2 t)) 

arises as a viscosity limit. We omit details, but remark that this can be shown 
along the lines of [I], where a similar analysis was carried out for the case of 
shear flows. 



6 Dissipative Solutions 

Let S(w) = |( Viu + Vu>') denote the symmetric gradient of a vectorfield w, and 



2 

set 



E(w) = -d t w - P(w ■ Vw), 



with P denoting the Leray-Helmholtz projection onto H(il). 

The following definition is from [52], given here in the version of [2] for 
bounded domains. The reader may consult these references also for a motivation 
of the definition. 

Definition 9. Let be a bounded domain with C 1 boundary. A vectorfield 
v € C([0,T];H w (Q)) is said to be a dissipative solution of the Euler equations 
^ if for every divergence-free test vectorfield w e C 1 (Q x [0, T] ) with wv \qq.= 
one has 

/ \v - w \ 2 dx<cxp\2 \\S(w)\\ tx> ds\ \v(x,0) -w(x,0)\ 2 dx 
Jn \ Jo J Jn ^25^ 

+2^ j cx p| 2 ^ \\S{w)\\ 00 dT^E{w)-{v-w)dxds 

for allie [0,T]. 

An immediate consequence of this definition is the weak-strong uniqueness 
(Proposition 4.1 in [22]): 

Proposition 10. Suppose there exists a solution v e C 1 (f2 x [0, T]) of the Euler 
equations (0). Then v is unique in the class of dissipative solutions with the 
same initial data. 

This follows simply by choosing w = v as a test function in the definition of 
dissipative solutions. 

Next, we prove Corollary [3] showing that admissible solutions may fail to be 
unique in bounded domains even for smooth initial data. 

Proof. Recall the construction from Section @] and define 

fi = {x € R 2 : p< \x\ < r Q } c fi. 
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It follows immediately from the definition that the restriction of a subsolution to 
a subdomain is itself a subsolution. Therefore we may consider the subsolution 
{v,u,q) constructed in Section U as a subsolution on f2 with energy density e as 
in Proposition with initial data given by 

x 1 

v(x,0) = -— r for x e ft 

(c.f. ((T210 . Applying this time Theorem [6] in O with this subsolution yields 
infinitely many admissible weak solutions as in the proof of Theorem [T] 

Since the initial data v(x,0) is smooth on f2, there exists a unique strong 
solution (indeed, this is the stationary solution). Thus weak-strong uniqueness 
fails, a fortiori implying that the non-stationary weak admissible solutions are 
not dissipative in the sense of Lions. □ 



7 A Criterion for Admissible Solutions to be 
Dissipative 

We have seen that, on bounded domains, an admissible weak solution may fail 
to be dissipative. However this will not happen provided such a solution is 
Holder continuous near the boundary of the domain, as claimed in Theorem 
[4] above. The aim of this last section is to prove this theorem. We follow 
Appendix B of [H] , but have to take into account that we need to deal with test 
functions which are not necessarily compactly supported in f2 in the definition 
of dissipative solutions. 

So let f2 be a bounded domain in M 2 with C 2 boundary and v an admissible 
weak solution of the Eulcr equations (|3|) as in the statement of Theorem [4] As- 
sume for the moment that for every divergence-free w e C 1 (fix [0,T]) satisfying 
the boundary condition we have 



dt 



v-wdx = / (S(w)(v - w) ■ (v - w) - E(w) ■ v) dx (26) 



in the sense of distributions, where E(w) is the quantity defined at the begin- 
ning of Section [51 We claim that (|26p implies already that v is a dissipative 
solution. Indeed this can be shown exactly as in : On the one hand, since v 
is admissible, 



I / 



in the sense of distributions. On the other hand, using the definition of E(w) 
and the identity J n (w ■ Vw) • wdx = (which follows from w ■ v Iqq= 0), we have 



dt 

Since 



,|2. 



dx = -2 / E(w) ■ wdx. (28) 



<> 



\v - w\ dx = / \v | dx + \w\ dx - 2 v ■ wdx, 
n Jn Jn Jn 
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we infer from this together with ([2"T|). and ((211) that 

— [ \v-w\ 2 dx<2 [ (E(w) ■ (v - w) - S(w)(v - w) ■ (v - w)) dx 
at Jn Jn 

<2 [ E(w) ■(v-w)dx + 2\\S(w)\\ oo [ \v - w\ 2 dx 
Jn °° Jn 

in the sense of distributions. We can then apply Gronwall's inequality as in 

to obtain (|25|) for every t € [0,T]. Therefore, it remains to prove ([26]) for every 

test function w. 

In [11] , identity ([26)) is proved for the case that w is compactly supported in fl 
at almost every time (see the considerations after equality (96) in fTT]). Let now 
w e C 1 (J1 x [0,T]) be a divergence-free vectorfield with w ■ = 0, which does 
not necessarily have compact support in space. We will suitably approximate w 
by vectorfields that do have compact support, much in the spirit of T. Kato [2T] 
(in particular Section 4 therein). 

Assume for the moment that f2 is simply connected, so that 90 has only 
one connected component. Since w is divergence- free, there exists a function 
ipe C([0,T];C* 2 (O)) nC'ffix [0,T]) such that 

w(x,t) = V L ip(x,t) 

and tp \dn= 0. Let now x '■ [0, °°) K be a nonnegative smooth function such 
that 

f if s < 1 
1 if s > 2 



X(s) 

and set 



Then, by Lemma 14.16 in [18], there exists rj > depending on f2 such that 
a; m- dist(x,90) is C 2 on 

T 7) = {xeTi:dist(x,dn) <T]}, (29) 

and hence w e e C*(f2 x [0,T]) for sufficiently small e > 0. Therefore, (|26p is true 
for w e : 

— f v-w e dx= [ (S(w e )(v - w e ) ■ (v — w e ) - E(w e ) -v) dx. (30) 
dt Jn Jn 

We will now let e tend to zero in order to recover (|2l)|) . 

Writing d(x) = dist(x,90), we have from the definition of w e : 



w e =X 



(^)^ + - e X'(~)^d, (31) 



and since ip e C([0,T];C 2 (O)) and tp \dn= 0, there is a constant C independent 
of t and e such that 

\if>(x,t)\<Cd(x) 
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for all x € Q,. Moreover, as the support of %'(-) is contained in (e,2e), and as 
|Vd| < 1, it follows from (J3TJ that 

w e ^w strongly in L°°([0,T];L 2 (CL)) (32) 

as e -» 0. For the left hand side of (|30|) this immediately implies 

— f v ■ w e dx -> — f v ■ wdx 
dt J n dt J n 

in the sense of distributions. Moreover, the right hand side of (f3"U| can be 
written, recalling the definition of E(w e ), as 

(S(w e )(v - w e ) ■ (v - w e ) - E(w e ) ■ v) dx 

[dfW f • v + (v ■ Vu> c ) • v - ((v - w e ) • Vw f ) • w e ] dx, 



and the right hand side of (|26[) is given by a similar expression. 
Next, observe that, again by (f3"2")l . 

dtw e ■ vdx -»■ / dtw ■ vdx 

in the sense of distributions and also that 

((v - w € ) ■ S7w e ) • w e dx = / ((v - w) ■ Vw) • wdx = 



thanks to the formula (v - w) ■ Vw) ■ w = (v - w) ■ ^V\w\ 2 and the fact that 
v - w e H(tt) (and similarly for ((u - w e ) ■ Vw e ) • w e ). 

To complete the proof of (|2l))) and therefore of Theorem^ it remains to show 
that 

/ (v ■ S7w e ) • vdx -*■ / (v-Vw)-vdx (33) 
Jn Jn 

in the sense of distributions as e -> 0. 

To this end, note that for every x e fl sufficiently close to <90 there exists a 

unique closest point x e d£l, and then 

x = x + d(x)v(x). 

We denote by t(x) = (-^(x), V\{x)) the unit vector at x tangent to d£l and use 
the notation v T (x) = v(x) -t(x), d T w v {x) = Vw v {x) -t{x), etc. (recall that x is 
uniquely determined by x). If e is sufficiently small, we can then write (recall 

m) 

(v ■ V(w t - w)) ■ vdx = / v v d v (w e - w) v v v dx + / v u d v {w t - w) T v T dx 
•/r 2e Jr 2 e 

v T d T (w e - w) v v u dx + I v T d T (w e - w) T v T dx 
Jr 2e 

-■ h +I 2 +h + Ii. 
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Recalling (f3~Tj) as well as V L ip = w an d observing Vd = v, we compute 

.(*(!)-!)*, 
(u> e - iy) T = - (x (-) " i) <^V> + -x' (-) V> 

^(lUe-uj)^ = -x'(-) w„ + (x(-) - iW^, (34) 
9„K - w) r = (x(~) - l) ^ + ^x" (^) V>, (35) 
9 T («; e -w)„ = (x(^) " ij <9 r w„, (36) 

9r(We -W) T = (x(^J - ij 9rW T + ^X'(^) w v (37) 

Before we estimate I\-Ia using (|34l) - (f37|) . let us collect some more information: 
As mentioned above, there is a constant C independent of t such that |^(a;)| < 
Cd(x). Moreover, since w e C 1 (r2 x [0,T]) and w v = on <9f2, we find similarly 
a constant independent of t such that |z/;„(:r)| < Cd{x). By assumption, if e is 
small enough, then v is Holder continuous with exponent a on r2 £ , uniformly 
in t, and since v e H(£l) implies v v = on <9fJ (cf. [16] Chapter III), we obtain 
another time- independent constant such that |i;„(x)| < Cd(x) a on I^. Finally 
note that ip, w, and v are uniformly bounded on I^ provided e is small, and 
that there is a constant independent of e such that |r2 £ | < Ce. 

In the light of these considerations we can use (|34[) - (j3"T)l to estimate 

\h\<- I vl\\ X '\L\\w T \\ 00 dx+ I «2|x-lLllft,uv|Lda; 
e Jr 2e Jr 2c 

< Ce 2a+1 + Ce 2a+1 -> 0, 

\h\< f MM ||x-l||oo \\d v w T \\ BB dx + — I KIKIIIx"lloolV'l^ 
< Ce a+1 + Ce Q -+ 0, 

\h\< [ Klklllx-i|LII<9r^L^<Ce Q+1 -o, 

Jr 2e 

\h\< [ v 2 T \\ x -i\L\\drW T \\ 00 dx + - [ ^Ix'ILKI 

< Ce + Ce ^ 0, 
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all estimates being uniform in time. This proves Theorem 0] if f2 is simply 
connected. 

As a final step, we convince ourselves that the proof can easily be modified to 
the general case when dfl has N connected components T 1 , . . . , T N in the spirit 
of Section 1.4 of [2D]. There still exists ip e C([0, T]; C 2 (f2)) n C 1 ^ x [0,T]) 
with V 1 ^ = w, but we can no longer require ip fao= 0. Instead, ip will take the 
constant value ip 1 on V, but the numbers ip 1 may be different. Now, if e > is 
small enough, then the sets 



with x as m the simply connected case. With this choice of w e we can then 



Remark 11. Theorem 0] implies that there can not be wild solutions on an 
annulus with smooth rotational initial data that are Holder continuous. Indeed, 
any admissible Holder continuous solution must be dissipative by our theorem, 
and the weak-strong uniqueness then yields that this solution must coincide with 
the stationary one. This observation is particularly interesting in the light of 
recent results (e.g. [DUDIITD] ) where examples of Holder continuous wild solutions 
are constructed. 

One of the first papers of Professor Mark Vishik "On general boundary 
problems for elliptic differential equations" \32j was essential, in particular in 
France, for the training of mathematicians in the generation of the first author 
of this contribution. Then when he turned to Navier-Stokes and turbulence he 
took an important role in progress over the last 60 years toward the mathemat- 
ical understanding of turbulence in fluid mechanics. Hence we hope that this 
essay will contribute to his memory and to the recognition of his influence on 
our community. 
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will be mutually disjoint, so that w c is well-defined by setting 




employ the very same arguments as above. 



□ 
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